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We study the phase structure of the four-dimensional twisted Eguchi-Kawai model using numer- 
ical simulations. This model is an effective tool for studying SU{N) gauge theory in the large-A^ 
limit and provides a nonperturbative formulation of the gauge theory on noncommutative spaces. 
Recently it was found that its symmetry, which is crucial for the validity of this model, can 
break spontaneously in the intermediate coupling region. We investigate in detail the symmetry 
breaking point from the weak coupling side. Our simulation results show that the continuum limit 
of this model cannot be taken. 
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1. Introduction 

In 1982, Eguchi and Kawai introduced an important and interesting idea, which is now called 
Eguchi-Kawai equivalence [jl|]. Consider the SU {N) gauge theory (YM) on aperiodic D-dimensional 
lattice with the Wilson plaquette action. In the large-A^ limit the space-time degrees of freedom 
can be neglected, and the theory is then equivalent to a model defined on a single hyper-cube, 
called the Eguchi-Kawai model (EK model). This correspondence was shown by observing that 
the Schwinger-Dyson equations for Wilson loops (loop equations) in both theories are the same. 
Naively, in the EK model the loop equations can have open Wilson lines, which do not exist in 
the original gauge theory due to gauge invariance. Therefore we need to assume that the global 
symmetry — e'^'' , which prohibits non-zero expectation values of the open Wilson lines, 
is not broken spontaneously. However, soon after the discovery of the equivalence, it was found 
that the symmetry is actually broken for Z) > 2 in the weak coupling region [Q]. Although the 
naive EK equivalence does not hold, several modifications have been proposed for this issue; they 
are the "quenched" Eguchi-Kawai model (QEK model) |Q ^, ^ and the "twisted" Eguchi-Kawai 
model (TEK model) [^. Historically, more work has been performed on the TEK model because 
it is theoretically interesting and numerically more practical. In addition, this model also describes 
gauge theories on noncommutative spaces (NCYM) |§, 0]. 

The TEK model is a matrix model defined by the partition function 

r ^ 

Ztek = Y\dU^, qx^{-Stek) (1.1) 

with the action 

Stek = -^N ^ Z^vTr U^U^Upl j8 = l/g", (1.2) 

where dU^ and {jJ. = I, are Haar measure and link variables. The phase factors Z^v are 

Z^v = exp {inin^v/N) , n^y = —ny^i G '^N- (1-3) 
The classical solution u'^^ = satisfies the 't Hooft algebra 

r^Ty = ZvjurvF^, (1-4) 

and is called "twist-eater". The twist-eater guarantees existence of the symmetry in the weak 
coupling limit. It is unclear whether or not the symmetry is unbroken in the intermediate coupling 
region, as there is no guarantee the symmetry is preserved. Numerical simulations performed 
during the early days of this model, however, showed that the symmetry is unbroken throughout 
the whole coupling region. This has encouraged the belief that the TEK model describes the large- 
N Umit of SU {N) gauge theory. 

Surprisingly, some indications of Z^ symmetry breaking were recently reported in several 
contexts concerning the TEK model [|8[ ^ |l^]. In [[lO|], the D = 4 TEK model with standard twist 
was investigated uptoN = 144 and Z^ symmetry breaking phenomena in the intermediate coupling 
region was observed by Monte-Carlo simulations. The work presented in this article continues this 
investigation. We concentrate on investigating the locations of the symmetry breaking from the 
weak coupling side in the {P,N) plane to determine if the continuum limits of the TEK models can 
be taken as the YM and the NCYM. 
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2. Twist prescriptions 

In this study we treat the D = 4 case. Among the various types of twist possible, we apply: 







A^ = 


L2 


(minimal symmetric twist, standard twist) 


liv 


T pSkew 


A^ = 




(minimal skew-diagonal twist), 


liv 


IflL Sjjy , 


N = 


ml} 


(generic skew-diagonal twist), 



(2.1) 
(2.2) 
(2.3) 



where we define anti-symmetrization matrices as 

/ 1 1 l\ 
-10 11 
-1-10 1 
V-1 -1 -10/ 



sym 



pSkew 



( ^ ^ 


o\ 


-1 








1 





-10^ 



(2.4) 



These twists represent lattices. The symmetric and the skew-diagonal form can be transformed 
into one another by an SL{A,'L) transformation for the coordinates on [11]. While these forms 



differ only by a coordinate transformation, they can give different results except the weak coupling 

ige theory on m-coincide 
1) of the generic twist.) 



hmit. We note that the generic twist (2.3) can be regarded as the gauge theory on m-coincident 
fuzzy T"^. (The minimal twists (2.1) and (2.2) are particular cases (m ■ 



3. Theoretical discussion for the symmetry breaking in the TEK model 

As we mentioned in the introduction, the symmetry can be broken in the intermediate 
coupling region. In this section we give a theoretical discussion about this phenomena. 

Here, we consider the first breaking point from the weak coupling limit jS^^. We assume that 
symmetry breaking at this point is a transition from the twist-eater phase = to the identity 

configuration phase = 1^} For simplicity we consider a type transition here. Of 



course we can treat Z^ Z^ — 



72 

-'N ' 



— ' Z^ (cascade) type breaking, but the obtained 
behavior is not different from the former type. First, we focus on the classical energy difference 



between these configurations. The energy difference can be easily calculated from the action ( |1.2D 
as 



^s = StekP^, = 1n]-Stek[u^ = r^] 

I \ A 



(3.1) 



For the generic twist we have 



AS: 



(symmetric form), 



(skew-diagonal form). 

'of course, the twist-eater only has symmetry for the twists we apply. But we write it as in this article. 



(3.2) 
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Note that the symmetric form is roughly three times more stable than the skew-diagonal form if both 
twists have similar quantum fluctuations. Thus the symmetry breaking for the skew-diagonal 
form can occur at smaller N than that for the symmetric form. 

Going away from the weak coupling limit, the system experiences greater quantum fluctu- 
ations. Here, we naively expect that the symmetry is broken if the fluctuation around the 
twist-eater configuration exceeds the energy difference AS'. Because the system describes 0{N'^) 
interacting gluons, it is natural to assume that their quantum fluctuations provide an 0{N^) con- 
tribution to the effective action. For the generic twist, the fluctuation is 0{m^L'^). Combining this 



with eq. (3.2), we can estimate the critical point as 



3L. .2 



(3.3) 



Although the above discussion is quite crude, the symmetry breaking behavior described by 



(3.3) is consistent with the numerical results discussed in the next section. 



4. Numerical simulations 

In this section we show the results of our numerical simulations for the symmetry breaking 
phenomena. In order to discuss the continuum and \wge.-N limits for this model, we concentrate 
on the first breaking point from the weak coupling side. 



4.1 Simulation method 



In our simulation we use the pseudo-heatbath algorithm. The algorithm is based on []12[], and 
in each sweep over-relaxation is performed five times after multiplying SU (2) subgroup matrices. 
The number of sweeps is 0(1000) for each j8. We scanned for the symmetry breaking point with 
a resolution of Aj3 = 0.005, and thus quote ±0.0025 as the error due to the finite resolution. Note 
that the breaking points are ambiguous because the breakdown of the symmetry is a first- 
order transition. As an order parameter for detecting the symmetry breakdown, we measure the 
expectation value of Polyakov lines 



1 

—TvU, 

N ' 



(4.1) 



4.2 Simulation results 



First we consider the minimal twists. Figures |I] and ^ show the critical lattice coupling from the 
weak coupUng side jS^ for the symmetric and skew-diagonal twists, respectively. For the minimal 
skew-diagonal twist we also observe the critical lattice coupling from strong coupling side j8^. 
We see that while the symmetry begins to break at N = 100 for the symmetric form, it is 
already violated at = 25 for the skew-diagonal form, which is consistent with the theoretical 
considerations in section ^. Additionally, we observe a clear linear dependence of jS/" on A'^(= L^): 

jS^^ ~ 0.001 1A^ + 0.21 (minimal symmetric twist), (4.2) 
j8c^ ~ 0.0034A^ + 0.25 (minimal skew-diagonal twist), (4.3) 
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o breaking point from vteak coupling side 
X breal<ing point from strong coupling side 




breaking region 



~2.9/N + 0.18 



50 



100 

N ( = l' ) 



150 



200 



Figure 1: j3^^ versus (minimal symmetric twist). Figure 2: fi^ and j3/^ versus (minimal skew- 
Fit line is eq. (42), which is obtained using > diagonal twist). Fit lines are eqs. (O) and {XA), 
169 data. which are obtained using > 64 and N>25 data, 

respectively. The symmetry is broken within 
the light blue shaded area. 



in the larger N region. This behavior was already obtained in the theoretical discussion. Note that 
the coefficient of N for the skew-diagonal twist is roughly three times larger than that for symmetric 
twist, which is also consistent with the theoretical analysis. For j8/^, we find the relation 

~ 2.9/A^ + 0.18 (minimal skew-diagonal twist). (4.4) 

As N is increased, approaches a point where the phase transition takes place in the 

original EK model. 

For the generic twist we use the skew-diagonal form because symmetry breaking occurs at 
smaller N than for the symmetric form, which makes our investigation much easier. We measure 
for this twist up to m = 4. The simulation results are plotted in figure |3[ From this figure we find 
that for each L, the show weak l/m behavior. The points at 1/m = are linearly extrapolated 
values, (m = oo means that an infinite number of fuzzy tori are superimposed.) The behavior for 
L = 5 is particularly interesting. While Zjj symmetry breaking is observed for m= 1,2, and 3, it is 
not seen for m = 4 because j3/' reaches a bulk transition point as m is increased. Figure ^ represents 
the same data as figure ||, but with as the horizontal axis. As we have seen in the minimal 
case, the data for L > 8 can be fitted by a linear function in for each m. From these figures, we 
discover that the data for L > 8 are well fitted globally by: 

„f , 0.060 

~ 0.0034L^ H h 0. 19 (generic skew-diagonal twist). (4.5) 

m 

The discussion in section ^ did not predict the observed dependence of on l/m. While we do 
not have a clear reason for this phenomenon at present, we suspect that it is related to collective 
modes of the eigenvalues of the link variables. 
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Figure 3: ^j; versus 1/m for each L (generic skew- Figure 4: )3,^ versus for each m (generic skew- 
diagonal twist). for m = oo is obtained by linear diagonal twist). We also include m = oo data, which 
extrapolation. is obtained by the extrapolation shown in figure ^ 



5. Continuum limit 



Although our simulation is restricted to the small N region, we may use our theoretical con- 
siderations to make statements about the large-A'^ limit. Thus the EK equivalence is valid only in 
the region jS > jS^ ~ A'^, not only in the smaller N region, but also in the \wge.-N limit. 

As is well known, the one-loop perturbative calculation of the YM lattice theory shows that its 
beta function behaves as j8 ~ logc?^^ near the weak coupling limit, where a is lattice spacing. If we 
wish to make the TEK model correspond to the YM theory, the scaling of the TEK model should 
obey that of the YM. In the TEK model, the lattice size L is related to N. (For the generic twist, the 
relation is A'^ = mL^.) Then the continuum limit of the YM system with fixed physical size I = aL 
can be obtained using the scaling 



j8 ~ log a ' ~ \ogN. 



(5.1) 



In order to obtain the \wge.-N limit with infinite volume, we should keep j8 lower than eq. (^). 
Otherwise, the system would shrink to a point. However, the simulation results obtained in this 
study show that grows faster than the logarithm. Therefore we conclude that EK equivalence 
breaks down and the TEK model does not have YM as its continuum limit. 

In the case of the NCYM, the beta function is essentially the same as that of the YM theory at 
the one-loop level [13|, and thus the scaling near the weak coupling limit is j8 ~ loga^^ But if we 
wish to make the TEK model correspond to the NCYM, there is the constraint a^L = fixed, which 
means that we take a scheme in which the noncommutative parameter Q is fixed. And then both 
the continuum limit and the infinite volume limit are simultaneously taken (double scaling limit). 
Regardless of the constraint, by the nature of the logarithm scaling, the scaling for the NCYM is 



the same as that of the ordinary YM (eq. 5.1). Therefore we conclude also that the TEK model 
does not have NCYM as its continuum limit. 
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6. Conclusions 

We carefully investigated the symmetry breaking phenomena in the TEK model using 
Monte-Carlo simulation. We found a clear linear dependence on for the symmetry breaking 
point from the weak coupling side. Regrettably, this means that the TEK models which use simple 
twists cannot be made to correspond to either ordinary YM or NCYM in the continuum limit. 

Finally, we mention the partial reduction p4]], which has been actively used in recent years. 



[14] showed that the \wge-N reduction is valid above some critical physical size Ic- This means 
that for a lattice size L the reduction holds only below some lattice coupling j8(L). In order to take 
continuum limit we should avoid the bulk transition at ^J? , and thus there is a lower limit to the 
lattice size L^. (j8) used for the reduction. It is clever that the twist prescription is applied to this 
reduction [lIS]], and we believe it would be very efficient. 



Acknowledgments 

The numerical computations for this work were carried out in part at the Yukawa Institute 
Computer Facility. M. H. is supported by Special Postdoctoral Researchers Program at RIKEN. 
T. H. would like to thank the Japan Society for the Promotion of Science for financial support. 

References 

[1] T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 1063. 

[2] G. Bhanot, U. M. Heller and H. Neuberger, Phys. Lett. B 113 (1982) 47. 

[3] G. Parisi, Phys. Lett. B 112 (1982) 463. 

[4] D. J. Gross and Y. Kitazawa, Nucl. Phys. B 206 (1982) 440. 

[5] A. Gonzalez-AiToyo and M. Okawa, Phys. Lett. B 120 (1983) 174; 

A. Gonzalez-Arroyo and M. Okawa, Phys. Rev. D 27 (1983) 2397. 
[6] H. Aoki, N. Ishibashi, S. Iso, H. Kawai, Y. Kitazawa and T. Tada, Nucl. Phys. B 565 (2000) 176. 

[7] J. Ambjorn, Y. M. Makeenko, J. Nishimura and R. J. Szabo, JHEP 9911 (1999) 029; 

J. Ambjorn, Y. M. Makeenko, J. Nishimura and R. J. Szabo, Phys. Lett. B 480 (2000) 399; 
J. Ambjorn, Y. M. Makeenko, J. Nishimura and R. J. Szabo, JHEP 0005 (2000) 023. 

[8] T. Ishikawa and M. Okawa, talk given by T. Ishikawa at the annual meeting of Japan Physics Society, 
Sendai, Japan, March 2003. 

[9] W. Bietenholz, J. Nishimura, Y. Susaki and J. Volkholz, JHEP 0610 (2006) 042. 

[10] M. Teper and H. Vairinhos, Phys. Lett. B 652 (2007) 359. 

[11] P van Baal and B. van Geemen, J. Math. Phys. 27 (1986) 455. 

[12] K. Fabricius and O. Haan, Phys. Lett. B 143 (1984) 459. 

[13] S. Minwalla, M. Van Raamsdonk and N. Seiberg, JHEP 0002 (2000) 020. 

[14] R. Narayanan and H. Neuberger, Phys. Rev. Lett. 91 (2003) 081601. 

[15] A. Gonzalez-Arroyo, R. Narayanan and H. Neuberger, Phys. Lett. B 631 (2005) 133. 



7 



